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ABSTRACT 

Theoretical  solutions  for  the  gap  thickness  in  an 
infinitely  wide  foil  bearing  are  reviewed.   An  experiment 
is  described  using  a  stainless  steel  foil  wrapped  about  a 
high  speed  rotor.   Details  of  instrumentation  are  included. 
Results  are  presented  for  the  relation  between  gap  thickness 
and  foil  stiffness  for  low  values  of  foil  tension. 

Only  a  small  amount  of  data  was  obtained  due  to  failure 
of  the  rotor  support  mechanism.   The  results  are  therefore 
inconclusive.   The  results  do,  however,  indicate  a  dependence 
of  foil  gap  thickness  on  foil  stiffness.   The  influence  of 
fluid  inertia  on  foil  gap  thickness  is  also  indicated 
although  this  effect  appears  to  be  small  compared  to  the 
effects  of  foil  stiffness. 
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1.   Introduction 

The  subject  of  this  study  is  a  foil  type  bearing.   Such 
a  bearing  consists  of  a  compliant  film  moving  relative  to  a 
rigid  cylinder.   Figure  1  illustrates  a  typical  foil  bearing 
schematically.   This  bearing  configuration  has  found  appli- 
cation in  the  high  speed  tape  transport  systems  of  digital 
computers,  in  the  manufacture  of  cloth  and  in  fact,  anywhere 
that  a  film  moves  over  a  surface  with  high  relative  speed. 

Of  particular  interest  in  a  study  of  the  foil  bearing 
is  the  thickness  of  the  lubricating  film;  that  is,  the  gap 
thickness  between  the  foil  and  the  cylinder.   Also  of  inter- 
est is  the  contour  of  the  foil  around  the  cylinder.   Intui- 
tively we  expect  that  the  gap  thickness  will  be  affected  by 
physical  conditions  characterized  by  the  following  physical 
quantities:   relative  speed  between  foil  and  cylinder;  foil 
tension  and  foil  stiffness;  and  lubricant  density,  viscosity, 
and  compressibility. 

Several  non-dimensional  parameters  have  been  devised  to 
indicate  the  relative  importance  of  specific  physical  effects 
on  the  foil  bearing.   The  foil  bearing  number  relates  the 
viscous  forces  in  the  lubricating  film  to  the  average  pressure 
above  atmospheric  in  the  film.   The  compressibility  parameter 
relates  the  ambient  pressure  to  the  average  pressure  above 
atmospheric  in  the  film.   The  inertia  parameter  relates  the 
fluid  inertia  forces  to  the  average  pressure  above  atmospheric 
in  the  film.   The  stiffness  parameter  relates  the  bending 
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moment  in  the  foil  to  the  moment  of  the  foil  tension  about 
the  cylinder  axis. 

A  mathematically  convenient  idealization  of  the  foil 
bearing  problem  can  be  obtained  by  making  the  following 
assumptions;  the  foil  is  infinitely  wide;  the  foil  is 
perfectly  flexible;  the  lubricating  gas  is  incompressible; 
inertia  forces  are  negligible.   These  assumptions  correspond 
to  setting  the  inverse  compressibility  parameter,  the 
stiffness  parameter  and  the  inertia  parameter  equal  to  zero. 
Several  analyses  have  been  made  with  these  assumptions  im- 
posed.  The  results  of  these  analyses  have  shoxm  that  the 
nominal  gap  width  is  equal  to  a  constant  times  the  foil 
bearing  number  to  the  2/3  power.   These  results  are  valid 
for  a  broad  range  of  foil  bearing  applications. 

The  purpose  of  this  experiment  is  to  investigate  the 
effects  of  stiffness  on  the  gap  width.   A  convenient  way  to 
include  the  effect  of  stiffness  is  to  consider  its  effect 
on  the  results  of  the  ideal  solution.   Since  a  theoretical 
study  of  this  nature  is  available,  it  can  be  used  for  com- 
parison with  experimental  results. 
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2.   Review  of  Literature 

A  comprehensive  review  of  the  literature  pertaining  to 
the  foil  bearing  is  presented  by  L.  Licht  (10).   Some  of 
the  articles  he  has  reviewed  are  mentioned  here  along  with 
some  others  that  have  been  published  recently. 

The  first  study  of  the  foil  type  bearing  was  an  experi- 
ment by  H.  Blok  and  J.  van  Rossum  (1).   They  also  included 
an  analytical  discussion.   They  present  a  solution  for  the 
nominal  gap  width  as  a  constant  times  the  foil  bearing  number 
to  the  2/3  power.   This  result  is  derived  from  the  Reynolds 
equation.   The  assumption  that  the  foil  contour  in  the  inlet 
region  is  parabolic  is  also  used.   Use  of  the  Reynolds 
equation  to  describe  the  fluid  behavior  assumes  negligible 
fluid  inertia  and  fluid  incompressibility. 

Cameron  and  Pat el  (2)  derived  a  less  restrictive  differ- 
ential equation  for  the  gap  thickness  than  Blok  and  van 
Rossum.   They  employ  the  Reynolds  equation  in  polar  coordi- 
nates.  The  pressure  term  in  this  equation  is  eliminated  by 
assuming  negligible  foil  stiffness.   In  this  case  equilibrium 
requires  the  fluid  pressure  to  be  equal  to  the  ratio  of  the 
foil  tension  to  the  local  radius  of  curvature.   Only  first 
order  terms  in  the  equation  for  the  radius  of  curvature  in 
polar  coordinates  are  retained.   They  did  not  present  a 
solution  to  this  equation. 

Baumeister  (3)  presents  a  solution  to  the  Cameron-patel 
equation.   This  solution  is  obtained  using  numerical  methods. 
He  also  presents  solutions  to  more  complicated  differential 
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equations  which  are  derived  in  the  manner  of  the  Cameron- 
Patel  equation  but  retain  second  and  higher  order  terms  in 
the  expression  for  the  local  radius  of  curvature.   He  does 
not  indicate  the  details  of  the  numerical  analysis. 

Eshel  and  Elrod  (5)  have  presented  a  more  refined  and 
elaborate  solution  for  the  nominal  gap  width  than  those 
already  mentioned.   They  employ  the  same  limiting  assumption 
as  Cameron  and  Pat el  and  Baumeister.   In  tnis  work  Eshel  and 
Elrod  utilize  a  small  perturbation  analysis  to  reduce  the 
Reynolds  equation  to  a  zeroth  order  differential  equation 
in  a  non-dimensional  film  thickness.   Details  of  the  numerical 
analysis  are  also  presented.   Extension  of  their  solution 
into  the  exit  region  indicates  the  presence  of  undulations 
in  the  foil  contour. 

Eshel  and  Elrod  (6)  have  studied  the  effects  of  foil 
stiffness  on  the  gap  width.   In  this  analysis  all  the 
assumptions  of  their  earlier  analysis  are  imposed  except 
that  of  negligible  foil  stiffness.   The  results  of  this 
analysis  are  presented  as  a  correction  curve  to  the  results 
obtained  in  their  more  restrictive  paper  (5). 

In  a  further  refinement  of  their  basic  study,  Eshel 
and  Elrod  (7)  have  studied  the  effects  of  compressibility 
on  the  nominal  gap  width.   Again  all  of  the  original  limiting 
assumptions  are  made  except  that  of  fluid  incompressibility. 
As  in  their  study  of  stiffness  effects,  the  results  are 
presented  as  a  correction  curve  to  the  results  presented  in 
the  first  paper  (5). 
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Several  experimental  studies  have  been  conducted  by 
L.  Licht  (8),  (9),  and  (10).   Licht  uses  a  small  capacitive 
probe  mounted  in  the  rotor  to  measure  the  gap  thickness 
directly.   Experiments  are  conducted  with  metalic  foils  as 
well  as  with  more  flexible  Mylar  foils  coated  with  a  con- 
ducting film.   Results  of  these  experiments  corroborate  the 
theoretical  results  of  Baumeister,  and  Eshel  and  Elrod. 
Licht  also  presents  oscillographs  which  show  the  undulations 
of  the  foil  in  the  exit  region. 

Experiments  made  by  Ma  (11)  also  are  in  agreement  with 
the  theory  of  Eshel  and  Elrod  and  Baumeister.   Ma  employed 
a  rectangular  probe  mounted  on  a  rotating  drum.   The  probe 
spanned  three-fourths  of  the  foil  width  and  measured  the 
average  gap  thickness  directly. 

Licht  and  Eshel  (16)  have  recently  presented  a  thorough 
study  of  a  foil  rotor  support  system.   A  high  speed  rotor  is 
supported  by  three  symetrically  located  foil  segments,  each 
of  itfhich  has  a  60  degree  angle  of  wrap  about  the  rotor. 
Foil  gap,  foil  contour,  and  response  of  the  system  to  various 
imposed  disturbances  are  studied.   Applicable  theoretical 
considerations  are  presented  in  detail. 

Eshel  (17)  has  recently  made  an  analytical  study  of 
the  influence  of  fluid  inertia  forces  on  the  nominal  gap 
width.   This  study  has  necessitated  a  new  reduction  of  the 
Navier-Stokes  equations.   In  this  analysis  he  has  also  in- 
cluded the  effects  of  fluid  compressibility.   Eshel  has 
reduced  the  Navier-Stokes  equations  to  prandtl's  boundary 
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layer  equations.   These  equations  are  solved  simultaneously- 
using  numerical  methods.   The  solution  is  presented  in  the 
form  of  correction  curves  to  the  results  obtained  in  his 
basic  study  (5). 
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3.   Theoretical  Considerations 

The  governing  equations  for  the  fluid  behavior  in  a 
foil  bearing,  as  in  many  fluid  problems,  are  the  Navier- 
Stokes  equations.   Several  physically  reasonable  assumptions 
pertaining  to  the  foil  bearing  can  be  made  that  will  reduce 
these  complicated  equations.   These  assumptions  are  listed 
here. 

(1)  The  foil  bearing,  is  infinitely  wide. 

(2)  The  fluid  flow  is  steady  and  one  dimensional. 

(3)  The  gap  thickness  is  small  compared  to  the  cylinder 
radius. 

(k)      The  fluid  is  incompressible. 

(5)   Inertia  forces  are  small  compared  to  viscous  forces. 
With  these  assumptions  imposed  the  governing  equation  of 
lubrication  for  the  foil  bearing  becomes  the  Reynolds 
equation. 

d/ds(h3dF/ds)  =  +  6  m  U  (dh/ds)  [l] 

h  =  the  thickness  of  the  gap  between  the 
foil  and  the  cylinder 

p  =  the  fluid  pressure 

s  =  the  distance  along  the  cylinder 

U  =  the  relative  velocity  between  the  foil 
and  the  cylinder 

M  -   the  fluid  viscosity 

Gross  (12)  presents  a  thorough  derivation  of  this 
equation  from  the  Navier-Stokes  equations. 
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Blok  and  van  Rossum  (1)  assumed  that 

h  =  h*  +  s2/2r  [2] 

where     h  =  the  nominal  gap  thickness  (where  dp/ds  =  0). 

They  obtained  the  following  equation  for  the  nominal  gap 
thickness. 

h*/rQ  =  .^26(6yau/T)2/3  [3] 

Cameron  and  patel  (2)  expressed  the  Reynolds  equation 
in  polar  coordinates. 

d/dQ(h3  dp/dO)   =  +  ^«.Ur(dh/dQ)  [>] 

h  -  r  -   rQ 

r  =  the  radial  coordinate 
0  =  the  angular  coordinate 

They  also  introduced  the  non-dimensional  coordinate 

X  =  h/rQ  . 

They  assumed  a  perfectly  flexible  foil  so  that  equilibrium 
of  the  foil  required  that 

P  -  PQ  =  T/R  [5] 

where     R  =  the  local  radius  of  curvature 

T  =  the  foil  tension  per  unit  width 

p  =  the  ambient  pressure 
o 
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The  radius  of  curvature  in  polar  coordinates  is 

H  =      [(dr/dG)2  +  r2]3/2  -. 

2(&r/&Q)z   -  r(d2r/d02)  +  r2  ' 

Expanding  this  relationship  for  curvature,  and  retaining  only 
terms  of  1st  order,  and  using  equation  [5]  to  eliminate  the 
pressure  term  from  the  Reynolds  equation,  a  non-dimensional 
differential  equation  is  obtained. 

d/dO(X  +  d2X/d02)  =  +  (6^U/T)(X  -  X*)/X3  [7 J 

X  =  the  nominal  non-dimensional  gap  thickness 

Baumeister  (3),  (M  has  solved  this  equation  using 
numerical  methods  and  has  obtained  the  following  equation. 

h*/rQ  =  .639  (6mU/T)2/3  [8 J 

He  also  has  solved  a  more  complicated  equation  derived  as 
above  but  retaining  second  and  higher  order  terms  in  the 
expansion  of  the  curvature  equation.   This  result  is  as 
follows : 

h*/rQ  =  .650  (6Mu/T)2/3  [9] 

Eshel  and  Elrod  (5)  have  provided  a  more  refined 
solution  for  the  nominal  gap  width.   The  assumptions  made  are 
the  same  as  those  mentioned  above.   The  derivation  of  the 
differential  equation  employs  a  similar  approach  to  that  of 
Cameron  and  patel  except  the  non-dimensional  film  thickness 
is  assumed  to  be 
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h/rQ  =  1(0  )en  [10 J 

where     11(6)  =  the  non-dimensional  gap  thickness 
6.  =  a  small  parameter 
n  =  an  undetermined  exponent 

h*/r0  =  H*en  . 
They  also  introduce  an  extended  coordinate 

I  =  e<rm  [ll] 

where     m  =  an  undetermined  exponent. 

In  addition,  a  non-dimensional  pressure  is  defined. 

tt  =  (*  -  P0)rQ/T  =  Tr0/TH  -  rQ/B  [12 J 

Using  these  equations  tt  can  be  expressed  as  a  series  in  €.  so 
that 

tt  «  1  +  (d2H/dJ2)€n-2m  +  (dH/d£)2/2  +  2H(dH/d^)e2n"2m 
-  H£n  +  3/2  (d2H/d£2  )  (dH/d£ )2£3n_2m  + [13J 

Equation  [4]  can  then  be  reduced  to 

d/dJ>(TI3dTr/dp£3n-2n  =  +    (6AU/T)    (dH/d^)6n"m  [14] 

By  applying   the   boundary  conditions   of   the   problem   to 
equations    [13 J   an<3-   [l^J»    they   show   that 

n  =   2/3  a  -   1/3 

£  =  6/4U/T,  the  foil  bearing  number  which  relates  viscous 

forces  to  the  average  pressure  above  ambient 

in  the  fluid. 
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H  is  then  expanded  into  a  series  in  €  ,  where 

H  =  Hq  +  H1€2/3  +  •  •  •  .  [15J 


h*  =  ^  +  H*e  2/3  + 


Using  equation  [13]  in  equation  [14]  to  eliminate  tt,  a  family 
of  differential  equations  is  obtained.   The  zeroth  order 
equation  is 

a3H0/a?3  =  +  (h0  -  i*)/al    .  [16] 

Since  €  is  a  small  parameter  it  is  reasonable  to  accept  the 
zeroth  order  approximation  for  H*.  Eshel  and  Elrod  obtain 
H*  =  .6^3  so  that  in  the  region  of  uniformity 

h*/rQ  =  ,6i^€2n   -  [17 J 

Eshel  and  Slrod  have  also  analyzed  the  effects  of 
stiffness  on  the  solution  for  the  gap  thickness.   The  follow- 
ing equations  must  be  introduced  for  the  balance  of  forces 
and  moments  of  a  foil  element.   These  equations  modify 
equation  [5]. 

dT/ds  +  Q/H  +  %    =   0  [18] 

P  _  P  =  -  aU2/H  +  T/R  -  dVds  [19] 

dM/ds  -  Q  =  0  [20] 

£  =  (-h/2r  )(dp/d0)  +   U/h,  the  fluid  shear 

traction  on  the  foil 

H  =  D/ii,  the  bending  moment  per  unit  width 

D  =  Et^/12(l  -V2),  the  flexural  rigidity  of  the 

foil  per  unit  width 
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E  =  the  modulus  of  elasticity 
t  =  the  foil  thickness 
y '=   the  poisson  ratio 

In  this  developement  the  stiffness  parameter  is  introduced 

This  parameter  relates  the  moment  in  the  foil  to  the  moment 
of  tie  tension  about  the  center  of  the  cylinder.   The  same 
method  of  analysis  is  employed  as  in  the  derivation  of  [17] 
above.   The  result  of  this  analysis  is  presented  as  a  cor- 
rection to  equation  [17]  in  the  form  of  curve,  Figure  2. 

In  ashel*s  most  recent  analysis  of  the  foil  bearing 
problem  (17),  he  includes  both  inertia  and  compressibility 
effects.   He  seeks  a  solution  to  the  Navier-Stokes  equations 
and  assumes  only  isothermal,  two-dimensional,  steady  flow 
and  a  perfectly  flexible  foil,  tie   introduces  tne  following 
non-dimensional  parameters. 

C  =  PQr0/T,  the  compressibility  parameter,  which  relates 
atmospheric  pressure  to  the  average  pressure 
above  atmospheric  under  the  foil 

I  =  nPU^r  /T,  the  inertia  parameter,  which  relates 


2V 

o' 

inertia  forces  to  the  average  pressure  above 


atmospheric  under  the  foil 

The  following  non-dimensional  variables  are  also  employed  to 
reduce  the  Navier-Stokes  equations. 
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v  =     Ur/u£"k 

u  =    ire/u 
r(=   (y/rQ)e-n 

Now  the  following  conditions  are  imposed  on  the  problem. 

(a)  The  expansion  of  equation  [6]  for  the  curvature  must 
retain  the  capability  of  satisfying  the  condition 
that  the  foil  becomes  flat  at  oo; 

(b)  Some  inertia  term  must  be  the  same  order  as  some 
viscous  term; 

(c)  Lateral  and  longitudinal  terms  in  the  continuity 
equation  will  be  of  the  same  order. 

jilth   these  assumptions  and  the  parameters  and  variables 
introduced  above,  Eshel  shows  that  the  Navier-Stokes  equation 
reduce  to  the  zeroth  order  of  approximation  to  prandtl's 
boundary  layer  equations. 

id  +  n/c)  u(jVa^)  +  v(d&£«  )  =  -4W»  +  1/12^  n/»i£ 

[21] 
^  ^(C  +  n)v]  -i-  \,   {(C  +  f)S]  =  0  [22] 


rr  -  1    ^2H/v?  2  [23J 


,Jhen  these  equations  are  solved  simultaneously  by  numerical 
methods,  a  family  of  curves  is  obtained  relating  H*,  I,  and 
C.   These  curves  are  corrections  for  the  solution  for  H  in 
equation  [17].   One  of  these  curves  is  presented  in  Figure  3* 
The  following  equation  has  been  fitted  to  those  curves  for 
the  region 
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H 


0   <1/C<1 

o  <  i<  .5 
=   .6^3  +   .2861  -i-  1.905I2  -   .183/C 


[24] 


24 


4.   Experimental  Procedure 

The  purpose  of  this  experiment  is  to  study  the  effect 
of  foil  stiffness  on  the  gap  thickness  in  an  air  lubricated, 
self-acting,  foil  bearing.   Eshel  (6)  has  shoxra.  that  the 
effects  of  stiffness  on  foil  gap  are  related  by  the  stiffness 
parameter,  s  =  D £ ~  '  3   /   Tr^  which  is  the  ratio  of  the  moment 
in  the  foil  to  the  moment  of  the  foil  tension  about  the  rotor 
center.   Then  for  fixed  tension,  the  stiffness  parameter 
increases  as  the  foil  bearing  number,  f  =  6/AU/T,  decreases. 
The  method  of  approach  is  to  fix  the  foil  tension  and  run 
an  air  driven,  cylindrical  rotor  up  to  high  speed.   As  the 
rotor  is  slowed,  the  stiffness  parameter  increases  due  to 
the  effect  of  the  rotor  speed  on  the  foil  bearing  number. 

The  rotor  is  two  inches  in  diameter  and  7.1  inches  long, 
Figure  7.   The  rotor  is  made  of  stainless  steel  and  the  4.2 
inch  long  center  section  and  end  faces  have  been  precision 
ground  to  a  number  four  finish.   Both  waviness  and  roughness 
are  critical.   One  inch  long,  pressurized  air  journal  bearings 
are  located  at  each  end  of  the  center  section  to  support  the 
rotor  and  fix  its  position.   pressurized  air  thrust  bearings 
are  located  at  each  end  face  of  the  rotor. 

A  two  inch  wide,  one  one-thousandth  inch  thick,  stainless 
steel  foil  is  wrapped  around  the  center  of  the  rotor,  passed 
over  guides,  and  weighted  to  provide  tensioning.   A  mechanism 
for  changing  the  angle  of  wrap  of  the  foil  is  provided.   An 
illustration  of  the  apparatus  is  provided  in  figure  4. 
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The  journal  bearings  are  mounted  in  large  7-3/^  inch 
diameter  plates.   A  Bently-Nevada  type  318  AJX  inductance 
probe  is  mounted  to  one  of  the  plates  and  faces  into  the 
rotor  halfway  between  the  two  plates.   The  plates  ride  on 
ball  bearings  which  are  mounted  to  the  support  base.   The 
foil  is  also  mounted  to  this  base.   This  arrangement  permits 
determination  of  the  foil  gap  width  as  well  as  scanning  of 
the  foil  contour. 

The  probe  measures  the  distance  to  the  foil.  Sien   the 
rotor  is  not  turning,  a  reference  reading  is  taken  with  the 
foil  in  contact  with  the  rotor.   Comparison  of  probe  readings 
with  the  reference  readings  for  any  velocity  gives  the  foil 
gap  thickness  subject  to  a  small  correction.   Reference 
readings  must  be  taken  for  each  foil  tension  because  a  change 
in  foil  tension  may  change  the  position  of  the  rotor  center- 
line. 

Additional  probes  are  mounted  at  either  end  of  the  rotor. 
The  output  from  these  probes  is  displayed  on  two  oscilloscopes, 
one  for  each  end  of  the  rotor.   These  displays  indicate  the 
orbit  of  the  rotor  centerline  and  any  translation  of  the 
center  from  the  reference  position. 

The  rotor  is  driven  by  air  turbines.   These  turbines  are 
symmetrically  located  at  either  end  of  the  rotor  (see  Figure 

7). 

Maximum  shaft  speed  is  about  800  rps  with  minimum  foil 
tension  of  .39  lb/in.   Pull  throttle  shaft  speed  for  a  tension 
of  7.5  lb/in  was  600  rps. 

Experimental  procedure  is  as  follows: 


26 


The  foil  is  tensioned  by  placing  appropriate  weights  in 
the  weight  pan.   The  reading  of  the  foil  probe  is  recorded. 
The  oscilloscope  traces  are  centered.   The  shaft  is  brought 
up  to  maximum  speed.   After  steady  state  is  attained  (con- 
stant shaft  speed)  the  foil  probe  reading  is  recorded.   The 
rotor  orbit  and  the  movement  of  the  rotor  center  is  recorded 
from  each  oscilloscope.   This  completes  the  readings  for  one 
speed.   The  rotor  is  now  slowed  to  a  new  speed  and  the  above 
procedure  is  repeated.   This  process  is  continued  until 
constant  shaft  speed  can  no  longer  be  maintained.   ;Jhen  this 
situation  occurs  the  data  run  is  complete  and  the  shaft  is 
stopped  and  the  reference  reading  re-checked.   A  new  weight 
is  then  added  to  the  weight  pan  and  the  above  procedure  is 
repeated. 


27 


5.   Instrumentation 

The  instrumentation  for  this  experiment  performs  three 
functions.  They  are  measurement  of  the  following:  (i)  foil 
gap  thickness;  (ii)  rotor  position;  and  (iii)  rotor  speed. 

Distance  measurements  are  made  with  Bently-Nevada  type 
318  AJX  probes.   These  probes  are  driven  by  Bently-Nevada 
model  3100  N  detector  drivers.   The  power  supply  for  the 
drivers  is  a  Systems  Research  Corporation  power  Supply  model 
3566. 

The  Bently-Nevada  probe-proximitor  (detector  driver) 
arrangement  provides  a  d-c  output  voltage  proportional  to  the 
distance  from  the  probe  to  a  conducting  surface.   The  distance- 
voltage  relationship  is  linear  for  output  voltages  from  il-V  to 
14V.   For  stainless  steel,  this  gives  a  range  of  about  .001 
to  .009  inches,  and  provides  resolution  to  the  nearest 
.000005  inches  with  a  sufficiently  accurate  voltage  sensing 
device.   For  details  on  probe  calibration,  see  Appendix  A. 

The  output  of  the  foil  gap  probe  was  displayed  on  a 
Hewlett  Packard  model  2010C  digital  voltmeter  capable  of 
measuring  to  the  nearest  .001  volts  at  99  volts  full  scale. 
This  voltmeter  integrates  voltage  over  time  periods  of  .01, 
.1  and  1  second.   The  .1  second  time  period  was  employed. 

Rotor  position  was  determined  by  four  other  probes 
mounted  90°  apart,  two  at  each  end  of  the  rotor.   The  output 
from  these  probes  was  displayed  on  two  Hewlett  Packard  model 
130  oscilloscopes.   The  probe  outputs  were  biased  with  7.5 
volt  d-c  batteries  so  that  a  zero  reading  could  be  attained 
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near  the  center  of  the  probe  linear  range. 

Saaft  speed  measurement  was  accomplished  using  a  Hewlett 
Packard  506a  optical  tachometer  pickup  driving  a  Hewlett 
Packard  model  522  electronic  counter.   This  system  is  capable 
of  measuring  from  5  rps  to  5»000  rps.   The  pickup  was 
accomplished  by  blackening  a  3.^  inch  i\ride  half -sect ion  of 
the  shaft  adjacent  to  the  turbine  buckets  at  one  end. 
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6.   Results 

The  results  of  the  experimental  work  are  presented  in 

tabular  form  in  Tables  1,  2,  3  and  4  and  in  graphical  form 

in  Figures  10,  11,  12,  pud   13.   Figure  10  is  a  plot  of  the 

gap  thickness  h  against  t.ie  peripheral  velocity  U  where 

the  foil  tension  T  has  been  held  fixed  for  each  run.   Figure 

11  is  a  plot  of  h  against  the  foil  bearing  number,  £  =  ___  , 

T 

to  the  2/3  power.   This  is  done  because  in  all  of  the  theo- 
retical works  of  Eshel  the  non-dimensional  gap  width  is 
assumed  to  be  a  power  series  in  this  parameter. 

Figure  12  is  a  plot  of  the  nominal  non-dimensional  gap 
thickness  H  against  the  stiffness  parameter  S.   Figure  13  is 
a  plot  of  H~  against  the  inertia  parameter  I.   This  curve 
was  constructed  because  the  inertia  parameter  varies  sig- 
nificantly throughout  the  range  of  data,  and  it  would  be 
useful  to  separate  the  effects  of  inertia  from  those  of 
stiffness. 

The  angle  of  wrap  of  the  foil  for  both  data  runs  was  60 
degrees.   Foil  tension  was  .39  lb/in  and  1.19  lb/in  for  the 
respective  data  runs.   Shaft  speed  ranged  from  500  rps  to 
100  rps. 


30 


7.   Discussion  of  Results 

The  results  of  this  experiment  are  incomplete  and  there- 
fore inconclusive.   Only  two  data  runs  were  completed  due  to 
failure  of  the  rotor  support  system. 

During  the  first  attempt  at  a  data  run,  the  foil 
partially  slipped  out  of  the  foil  bracket.   As  a  result,  the 
effective  width  of  the  foil  went  from  two  inches  to  less 
than  ^  inch.   This  destroyed  the  air  film  between  foil  and 
rotor  and  in  turn  caused  one  of  the  bearings  to  fail.   Due 
to  the  large  dissipation  of  energy  in  such  a  short  time, 
the  bearing  seized  badly.   It  could  not  be  pressed  off  and 
had  to  be  cut  off  on  the  lathe. 

The  rotor  was  re-ground  and  a  new  bearing  was  manu- 
factured.  The  bearing  could  not  be  ground  as  the  proper 
facilities  were  not  available.   The  bearing  was  finished  by 
using  a  fine  grade  of  emery  paper.   This  refinished  rotor  and 
bearing  were  assembled  and  the  rotor  was  found  to  be  out-of- 
round. 

The  best  alternative  at  this  point  seemed  to  be  to  make 
a  new  rotor  and  to  have  the  existing  bearings  re-ground. 
Re-grinding  the  rotor  was  deemed  impractical  because  it 
necessitated  centerless  grinding.   No  centerless  grinding 
facilities  were  available. 

The  data  was  obtained  using  this  second  rotor,   pre- 
liminary testing  of  this  rotor  and  support  system  Indicated 
considerable  imbalance  but  with  a  very  low  resonance  point. 
As  time  was  very  short  it  was  decided  to  press  on  and  obtain 
some  data. 
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The  first  two  runs  of  data  were  made  without  difficulty. 
A  third  run  of  data  was  attempted  but  a  higher  rotor  speed 
was  necessary.   This  is  true  because  as  the  tension  is 
increased,  the  maximum  speed  must  also  be  increased  to  main- 
tain the  same  range  of  foil  bearing  number,   £  =  6/iXU/T. 
Each  succeeding  run  used  a  higher  foil  tension  than  the  one 
before. 

On  the  third  run  synchronous  whirl  became  evident  at 
about  650  rps.   Although  the  whirl  orbit  was  increasing  in 
size  as  the  speed  increased  it  was  nevertheless  small  compared 
to  the  clearance  in  the  bearing.   At  about  840  rps  both 
journal  bearings  failed  suddenly.   It  appears  that  a  resonance 
occurred  at  or  near  this  sp.od.   The  rotor  was  being  acceler- 
ated slowly  and  although  tae  whirl  orbit  was  still  growing 
slowly  at  about  800  rps^  it  seemed  to  be  stabalizing. 

Both  runs  of  data  that  were  obtained  were  runs  made  with 
constant  values  of  foil  tension.   The  inverse  compressibility 
parameter,  1/C  =  T/P  r  ,  is  constant  for  both  these  cases. 
Therefore  it  is  reasonable  to  neglect  the  compressibility 
effect  as  it  pertains  to  one  set  of  data. 

A  wide  range  of  stiffness  parameters  is  spanned  by  both 
runs  of  data,  see  Tables  2  and  4.   Unfortunately  a  wide  range 
of  inertia  parameters  is  also  spanned.   Recall  that  for  fixed 
foil  tension,  the  inertia  parameter  varies  as  velocity  to 
the  second  power.   For  this  reason,  it  is  impossible  to 
separate  the  two  effects. 

The  value  of  H*  =  .6^3  is  well  accepted  in  the  foil 
bearing  literature.   The  effect  of  inertia  is  to  increase 
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this  value,  while  compressibility  and  stiffness  tend  to 
decrease  this  value.   Both  runs  of  data  indicate  that  the 
effect  of  stiffness  for  this  range  of  speed  and  tension  is 
large  compared  to  the  effect  of  inertia. 

Comparison  between  the  two  runs  indicates  the  effect 
of  compressibility  over  a  wide  range  of  stiffness.   This 
indication  is  that  the  effect  of  compressibility  is  signifi- 
cant only  for  small  values  of  stiffness. 

These  results  are  more  useful  in  determining  the  proper 
conditions  for  further  research  than  they  are  in  giving 
quantitative  inference  into  the  behavior  of  a  foil  bearing. 
They  show  the  reader  that  in  order  to  study  the  effects  of 
compressibility  or  inertia,  the  stiffness  effect  must  be 
minimized. 
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8.   Recommendations  for  Further  iork 

Due  to  the  fact  that  so  little  data  was  obtained,  there 
seems  to  be  considerable  room  for  further  work  with  this 
apparatus.   Several  changes  need  to  be  made  before  the  device 
will  work  reliably  and  before  more  useful  data  can  be  ob- 
tained. 

The  significant  stiffness  effects  have  already  been 
mentioned.   A  good  project  would  be  to  eliminate  these  effects 
and  attempt  to  isolate  the  effects  of  Inertia  or  compressi- 
bility.  A  more  flexible  foil  such  as  Mylar,  coated  with  a 
thin  conducting  film,  could  be  tried  for  this  purpose. 

Also  mentioned  previously,  the  rotor  support  system, 
specifically  the  journal  bearings,  gave  considerable  trouble. 
The  bearings  were  made  out  of  stainless  steel,  the  same 
material  as  the  rotor.   This  was  done  to  minimize  temperature 
effects.   Consequently,  when  the  bearings  failed,  consider- 
able damage  was  done  to  both  bearings  and  journals. 

A.  necessary  improvement  would  then  be  the  re-design 
of  the  bearings  so  as  to  provide  protection  for  the  rotor  in 
the  event  of  failure.   To  this  end,  a  material  must  be  found 
that  is  soft  compared  to  stainless  steel,  yet  can  be  machined 
to  the  smoothness  required  of  an  air  bearing.   One  possi- 
bility might  be  to  re-machine  the  existing  bearings  and  coat 
the  bearing  surface  with  a  thin  film  of  Teflon.   Another 
possibility  might  be  to  use  nylon  bearings. 

The  cause  of  the  last  bearing  failure  was  probably 
Imbalance  of  the  rotor.   It  is  therefore  recommended  that 
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the  rotor  be  precision  balanced.   The  rotor  run-out  should 
be  held  to  less  than  .0001  inches  at  speeds  up  to  60,000  rpm. 
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TABLE  1. 

Measured  Quantities,  1st  Run 

Foil               Shaft  Gap 

Tension            Speed  Thickness 

lb/in              rps  volts 

.39              452  .735 

'128  .700 

41 6  .670 

393  .623 

360  .550 

300  .428 

250  .319 

199  .230 

152  .140 

101  .058 
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inches 

X  10< 

.000735 

118 

.000700 

112 

.000670 

109 

.000625 

103 

.000550 

95 

.000428 

78 

.000319 

66 

.000230 

52 

.000140 

39 

.000058 

26 

TABLE  2 

Computed  Quantities,  1st  Run 

Gap     Foil  Bear-    Inertia   Stiffness     Inverse 
Thickness   ing  Number  Parameter  Parameter   Compressibility 

Parameter 

1.1  2.93         .026 

.98  3.04 

.93  3.08 

.82  3.20 

.69  3.41 

.48  3. 81 

.34  4.33 

.21  5.06 

.12  6.06 

.05  7.80 
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TABLE  3. 
Measured  Quantities,  2nd  Run 


Gap 

Thickness 

volts 

.730* 
.361 

.333 
.282 

.25^ 
.220 
.177 
.150 
.133 


*  Tnese  measurements  made  at  an  earlier  time  than  the  rest. 


Foil 
Tension 
lb/in 

Shaft 

Speed 

rps 

1.19  mm 

775* 

492 

451 

^Ck 

350 

300 

242 

189 

160 

'!( 


TABLE  4, 

1 

poi: 

ing 
X 

Computed  Quantltl< 

is ,  2nd  Run 

Gap 
Thickness 
Inches 

L  Bear- 
Number 
105 

Inertia 
Parameter 

Stiffness 
Parameter 

.000730* 

67* 

1.04* 

1.40* 

.000361 

42 

.42 

1.88 

.000333 

39 

.36 

2.00 

.000282 

35 

.28 

2.16 

.000254 

30 

.21 

2.39 

.000220 

26 

.16 

2.68 

.000177 

21 

.10 

3.03  . 

.000150 

15 

.07 

3.78 

.000133 

14 

.05 

4.01 

Inverse 
unpressibil 
Parameter 

.08 


*  These  measurements  made  at  an  earlier  time  than  the  rest. 
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Foil  Tension  T 


Local  Gap 
Thickness 


h* 

Nominal  Gap 

Thickness 


FIGURE  1.  SCHEMATIC  OF  FOIL  TYPE  BEARING 
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FIGURE  2.  NOMINAL  CLEARANCE  VERSUS  STIFFNESS  AND  COMPRESSIBILITY  PARAMETERS 
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FIGURE  3.  NOMINAL  CLEARANCE  VERSUS  INERTIA  PARAMETER 

(from  Eshel  (17)) 
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PROBE  CALIBRATION  APPENDIX 

The  probe  calibration  system  consisted  of  a  Johannson 
Mikrckator   model  S510E-7  comparator  which  is  accurate  to 
the  nearest  .000005  inch.   The  comparator  was  fitted  with  a 
micrometer  adjustment  mechanism  salvaged  from  a  microscope. 
A  fixture  for  holding  a  probe  was  attached  to  this  mechanism. 
The  calibration  disc  in  the  comparator  was  made  of  the  same 
material  as  the  rotor. 

Moving  the  probe  up  and  down  with  the  micrometer  adjust- 
ment and  recording  the  voltage  output  on  the  digital  voltmeter 
provided  a  calibration  curve  for  a  small  range  of  voltage. 
This  was  done  for  several  nominal  voltages  for  each  probe- 
driver  combination.   These  curves  establish  the  linear 
operating  range  of  the  probe  and  give  the  slope  of  the  voltage 
versus  distance  line  for  that  range. 

The  probes  were  checked  for  temperature  effects  and  the 
following  conclusions  have  been  reached:  (i)   The  slope  of 
the  voltage  versus  distance  curve  is  not  temperature 
sensitive.   (li)   The  voltage  reading  for  a  given  distance 
is  temperature  sensitive.   (iii)   The  drivers  are  temperature 
sensitive,  but  only  for  large  variations  of  temperature. 

These  effects  are  not  considered  to  effect  significantly 
the  results  since  the  probes  are  located  in  positions  where 
large  temperature  variations  do  not  occur.   Driver  (ambient) 
temperature  is  constant  for  a  given  data  run,  and  varies  onl^ 
slightly  (less  than  10°)  over  all  runs.  A  calibration  curve 
has  been  provided  in  Figure  Ik, 


Calibration  for  the  foil  gap  measurement  was  accom- 
plished in  the  folio  ling  manner.   Several  thin  paper  shims 
were  measured  for  t  lickness  with  one  of  the  probes  mounted 
in  the  comparator.   The  shims  were  then  inserted  between  the 
foil  and  t  l^  rotor.   The  reading  on  the  digital  voltmeter 

-s  then  recorded.   A.  shim  was  then  romoved  and  another 
i'.  ding  taken.   This  procedure  was  continued  until  the  gap 
measuring  probe  was  out  of  its  linear  range  of  operation. 
The  foils  were  then  inserted  again  in  turn  and  readings 
were  taken  until  the  original  number  were  in  place.   Prom 
this  data  a  curve  was  obtained  for  gap  thickness  versus 
probe  output  voltage  differential.   The  resulting  calibration 
curve  is  presented  as  .Figure  16. 
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